Quantum Description of Nuclear Spin Cooling in a Quantum Dot 
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We study theoretically the cooling of an ensemble of nuclear spins coupled to the spin of a 
localized electron in a quantum dot. We obtain a master equation for the state of the nuclear 
spins interacting with a sequence of polarized electrons that allows us to study quantitatively the 
cooling process including the effect of nuclear spin coherences, which can lead to "dark states" of 
the nuclear system in which further cooling is inhibited. We show that the inhomogeneous Knight 
field mitigates this effect strongly and that the remaining dark state limitations can be overcome by 
very few shifts of the electron wave function, allowing for cooling far beyond the dark state limit. 
Numerical integration of the master equation indicates, that polarizations larger than 90% can be 
achieved within a millisecond timescale. 

PACS numbers: 71.70.Jp, 73. 21. La 
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I. INTRODUCTION 

Nuclear spins are one of the best studied quantum 
systems and highly developed techniques such as NMR 
have allowed detailed study of properties and dynamics 
of molecular and solid state systems Due to their 
very long decoherence time nuclear spins (and hyperfine 
levels) have also played a central role in many approaches 
to the implementation of quantum information process- 
ing (QIP) 

Recently, the localized ensemble of nuclear spins in a 
quantum dot (QD) has received special attention in the 
context of QIP with electron spins in QDs: the nuclei 
couple via a Fermi contact interaction to the electron 
spin and, as predicted by theory H, H, [M [H [11, 
have been shown in recent experiments to constitute the 
major source of decoherence of electron spin qubits in 
some of the most promising QD-based implementations 
[isl . [T^ . The vice of this strong coupling is turned into a 
virtue when the electron is used to manipulate the state 
of the nuclear ensemble. This has long been exploited 
in dynamical nuclear polarization (DNP) [l^ [11, |T3, Ull 
in bulk systems and afforded many insights in the spin 
dynamics in solids [13, [H ■ 

DNP in quantum dots has come into focus more re- 
cently in the context of QIP, since strongly polarized nu- 
clei could lead to much longer electron spin dephasing 
times [ll, provide strong local magnetic field gradients 
required in quantum information proposals [20l |21||. and 
even allow to utilize the nuclear spins themselves as long- 
lived quantum memory [13, HI] . More generally, a highly 
polarized nuclear spin ensemble in a QD provides, to- 
gether with the electron spin, a strongly coupled, well 
isolated mesoscopic quantum systems with close similar- 
ities to the Jaynes-Cummings model in quantum optics 
[23I . [2^ . [25j , with the fully polarized state corresponding 
to the vacuum in all cavity modes. Thus ultra- high DNP 
in QDs may open the door to realize cavity-QED in quan- 
tum dots and implement tasks such as state engineering. 

Experimentally, significant nuclear polarization in self- 
assembled QDs has been achieved [gl, [H [H [H [ij. 



However, the degree of polarization in these experiments 
was still too low to improve electron spin coherence times 
considerably and still far from the ground state. 

Theoretically, cooling dynamics has mostly been con- 
sidered in the spin temperature approximation [l|, 

[13, 

[sol . [Blj , in which coherences among the nuclear spins are 
neglected. This is appropriate if, as in bulk or quantum 
well systems, there is no fixed electron wave function and 
many motional states are involved, or if the nuclear de- 
phasing rate is large. In quantum dots, however, the 
nuclei interact collectively with an electron in the mo- 
tional ground state of the QD and the higher motional 
levels are far detuned. Therefore the coupling strength 
of each nucleus is fixed, and well defined phase relation- 
ships between the nuclear spins can build up, necessitat- 
ing a quantum treatment of the process, which was first 
pointed out by Imamoglu et al. [32l | , who showed that the 
cooling process can be inhibited by so-called dark states, 
which trap excitations and potentially result in serious 
constraints on the achievable polarizations. While it was 
pointed out in [33 | that inhomogeneities (either inherent 
in the system or introduced actively by modulating the 
wave function of the electron) can mitigate this prob- 
lem, these ideas were put to numerical test only in very 
small ID systems of 10 nuclear spins. However, the effect 
of inhomogeneities is expected to be reduced for realis- 
tic larger systems [22j, and thus limitations due to dark 
states are more severe (6Ci| . 

We consider the cooling of N nuclear spins in a QD 
through interaction with polarized electrons. One cooling 
cycle consists of (a) initialization of the electron spin in a 
well-defined direction, and (6) evolution of the combined 
system for a "short" time. In this way the electron spin 
acts effectively as a T = -reservoir for the nuclear spin 
bath, and pumps excitation out of it. 

We derive in a consistent manner a full quantum model 
of this process, which allows us to numerically study par- 
ticle numbers of up to ~ 10^. We show that a sufficient 
inhomogeneity of the couplings leads to a dephasing of 
nuclear spin states and thus limitations due to dark states 
are partially lifted. We demonstrate that enhanced cool- 
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ing protocols involving only a few (< 10) modulations of 
the electron wave function, allow to fully overcome these 
limitations, indicating that Overhauser fields above 90% 
of the maximal value can be created within the nuclear 
spin diffusion time. 

The paper is organized as follows: In SecHIlwe present 
the generic cooling protocol and analyze its performance 
in Sec. llin the applicability of the scheme to some specific 
physical systems is studied in Sec. |lVl 



II. THE COOLING SCHEME 

Interaction- The Fermi contact interaction between 
an (s-type conduction band) electron spin S and the spins 
li of the lattice nuclei leads to a Heisenberg like coupling 
Aaili ■ S to the nuclear spin at lattice site i, where A sets 
the overall strength of the hyperfine interaction and the 
factor < < 1 is determined by the probability to find 
the electron at site i and the gyromagnetic ratio of the 
ith nucleus fi\. In the presence of an external magnetic 
field Bcxt we write the Hamiltonian of the spin system 
with the collective nuclear spin operators — gil^ 
(/i — ±, z) as (ft = 1) 



H = I {A+S- 



S+A-) + gA'S' + .g>BSe.t5^ (1) 



where we have defined g — Ay'"^^ af and gi — 

Oiil \J^i a^, such that J2i9i = ^^'^ denoted the elec- 
tron factor by g* and the Bohr magneton by /is. 

We do not consider the Zccman energy of the nuclear 
spins, because for typical QDs it is much (10^ times) 
smaller than the electron's Zeeman energy and simi- 
larly we neglect the even smaller dipolar interaction be- 
tween the nuclei. The effects of these are briefly discussed 
at the end of Sec. IIIII Finally, we restrict the analysis to 
nuclear spins 1=1/2 and one nuclear species only in this 
article. 

The first part of the above Hamiltonian exchanges spin 
excitation between the electron and the nuclei, and it is 
this mechanism that is used to create polarization. The 
second part of the Hamiltonian constitutes a "quantum" 
magnetic field, the Overhauser field, for the electron spin 
generated by the nuclei. 

The cooling scheme- We assume initially the electron 
spin to be pointing in the — z-direction \ipe-) = li)- In 
the absence of a magnetic field this initial state defines 
the axis of quantization. The cooling cycle we consider is 
an iteration between evolution with Hamiltonian Eq. ([T]) , 
and reinitialization of the electron to ||). The nuclei ef- 
fectively "see" a large cold reservoir of electron spins and 
the concatenated evolution of the nuclear spin density 
matrix becomes 



P- 



■ Utile 



Ut{p®\i){i\)u^t]®\i)ml-- ■ (2) 



the following p will denote the state of the nuclear spin 
system only. Spin polarized currents or optical pump- 
ing with polarized light give rise to a polarized electron 
bath, but also the fast electrical control available in dou- 
ble QDs JJJ allows for the creation of nuclear spin po- 
larization without the need for pre-prepared electrons, as 
we will detail in the last section of this article. 

Considering small times for the evolution in each indi- 
vidual step of the cooling protocol, we expand the time 
evolution operators in Eq.® to second order. The stan- 
dard deviation of the A^'^^^-terms scales as A^f^~o?^ = 
g ^ 0{A/^/N) for the initially totally mixed nuclear spin 
state, and thus for At g^^ ~ ^/N j A we neglect higher 
orders. The readily obtained master equation 

.gAt, ^ g^{Mf . , 



{A+A-pt + ptA+A- - 2A-ptA+) , (3) 



Here Ut = exp{—iHt) is the time evolution operator, 
tr. denotes the trace over the electron, and here and in 



contains a Hamiltonian part arising from the Overhauser 
field and a contribution in Lindblad form. The latter gen- 
erates the nuclear spin polarization, and has been stud- 
ied in the limit of homogeneous coupling constants in the 
context of superradiance [H, [H, [3a| . 

As polarization builds up and ^(A^) :$> A/^/N the 
Hamiltonian terms on the right hand side of Eq. ([3]) may 
become large (for fixed time step At). To preserve va- 
lidity of the master equation one can either reduce the 
interaction time At < A~-^ or assume that the Over- 
hauser field (A^) is approximately compensated by an 
applied magnetic field, so that {gA^ — g* p,BBext)At <C 1 
for all times. In the latter case At is short enough to en- 
sure quasi-resonant hyperfine flips despite the random 
detunings stemming from the fluctuating Overhauser 
field and at the same time large enough to guarantee 
a fast cooling rate [6l|. This is the situation we in- 
vestigate in the following. Without retuning the sys- 
tem in this manner the polarization rate becomes de- 
pendent on the polarization itself and the emerging non- 
linearities give rise to the bistability effects observed in 
^14, 29, 36, 37.jj18, 39, 48]and limit the final polarization. 

Homogeneous Coupling- Before we discuss general in- 
homogeneous couplings, consider for a moment the ho- 
mogeneous case, ai oc as a demonstration of some 
interesting features of the above master equation. In this 
case, the operators A±^z appearing in Eq. ([3]) form a spin 
algebra I±^z and the collective angular momentum states 
(Dicke states) | J, m/, 0) provide an efficient description 
of the system dynamics |40| : the total spin quantum 
number / is not changed by A±^z and the effect of Eq. ([3]) 
is simply to lower (at an (/, TO/)-dependent rate) the 
quantum number. If m/ = — / is reached, the system can 
not be cooled any further, even if (for / <C N/2) it is far 
from being fully polarized. These dark states [13, [13] are 
a consequence of the collective interaction Eq. ^ . Thus 
spin excitations are trapped and cooling to the ground 
state prevented. We evaluate the steady state polariza- 
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FIG. 1: (Color online) Exact polarization dynamics. Left: 
Homogeneous case, gj = 1/VTl. Right: In the inhomoge- 
neous case, gj oc exp{—{j — N/2— /w^). The term 1/4 is 
added to account for asymmetry between electron wave func- 
tion and the lattice and avoid symmetry effects for this small 
scale system. 



tion (/^),s-(E.^f/ViV)- as 

V^-2«]vi:/P^ + l)'''-Vi^+°WiV). (4) 

i.e. for a mesoscopic number of particles the obtained 
polarization is negligible. In the above equation (/^)o is 
the expectation value in the completely polarized state, 

Di^(^ ^ I -l) '^^^^^ degeneracy 

of the subspaces of different total angular momentum, 
and the last equality has been obtained by employing 
the Stirling formula. 

Evolving the nuclei according to Eq.([3]), we find the ex- 
act time evolution of the polarization as shown in Fig. [T] 
In these and the following simulations gAt — 0.1, i.e. 
At = 0.1g~^ ~ Q.l\fN/A. As expected the polariza- 
tion decreases as as N increases, which underlines 
the importance of the nuclear spin coherences. In par- 
ticular this shows that an incoherent spin temperature 
description of the process would give even qualitatively 
wrong results. The timescale over which the steady state 
is reached is ~ N / {g/SiA). 

Inhomogeneous Coupling- Consider now an inhomoge- 
neous wave function. The results for the exact evolution 
of the quantity of interest, (A^), are shown in Fig.[T] The 
coupling constants gj in this example are taken from a 
ID Gaussian distribution with width iV/4. [62] The most 
important and striking feature is that in this situation 
almost complete polarization is obtained. 

The reason that this is possible here is not that there 
are no dark states in the case of inhomogeneous cou- 
pling constants. On the contrary it has been shown 
that there exists an one-to-one mapping [2^ from the 
famihar homogeneous dark states {|/, — /, (3) in the Dicke 
basis) to their inhomogeneous counterparts, defined by 
A^\D) = 0. The reason for obtaining high polarization 
beyond the homogeneous limit is the Hamiltonian part 
of the master equation ([3]). To illustrate this point, con- 



sider two spins with coupling constants gi ^ g2- Then 
the dark state \^ d) oc g2|Ti) — 5i|iT) evolves due to the 
A^-term in Eq.© to e**<'*g2|Ti)-e-*'^'*.gi|it), where Sg is 
proportional to (?i — (?2- Obviously this state will become 
"bright" again after a time cx l/\gi — gj \ and A^\D) ^ 0. 
This process is first order and, as we will detail later, "de- 
livers" coolable excitations sufhciently fast to maintain a 
high coohng rate. 

III. POLARIZATION DYNAMICS 

The polarization dynamics of the nuclear ensemble is 
governed by Eq. ([3]). While for homogeneous systems 
the collective angular momentum Dicke basis enables an 
efficient description of the problem, for realistic large and 
inhomogeneous systems more effort is required. 

To study the evolution of the nuclear polarization, we 
are interested in the individual spin expectation values 
(ct^cTj^). These depend, via Eq. ^ on all the elements 
of the covariance matrix 

which, in turn, depend on higher order correlations as 
seen from the equations of motion 

k 

+g,(a+[<yJ,cT+]a^)), (5) 

where = ig{gj — gi)/2 — g'^At{gj — gi)^/8 and k = 
g^At/8 and the erf refer to the Pauli matrices at site i. 

The simultaneous solution of the ensuing hierarchy of 
equations is only feasible for very small particle numbers 
N and further approximations are needed to treat the 
large systems of interest. We introduce several ways, 
labeled (i) to (v), of closing this set of equations and 
discuss their validity and implications in detail below. 

In the strongest approximation (i) all coherences be- 
tween different spins are neglected yielding independent 
rate equations for each individual nuclear spin. This 
reproduces essentially the spin-temperature description 
commonly employed in the discussion of bulk DNP P, [l3| 
(each subset of spins with identical coupling strengths gi 
is assigned its own effective temperature). This approach 
cannot reproduce the quantum effects we want to study, 
but it can serve as a benchmark for how strongly these 
are influencing the cooling process. 

The simplest approximations that take quantum co- 
herences between nuclear spins into account close the 
hierarchy of equations at the level of second order cor- 
relations. Our approximation (ii) is motivated by the 
generalized Holstein-Primakoff description ^4l]| , which in 
lowest order treats the nuclei as bosonic modes a~ — > a^. 
The bosonic commutations relations [ai,aj] = 6ij yield 
a closed set of equations for the elements of the covari- 
ance matrix 7. The bosonic description is known to be 
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FIG. 2; (Color online) Comparison of different approximation 
schemes for the homogeneous situation with A'' = 100 (left) 
and the case of Gaussian couplings (as in Fig. [1]) and = 10 
nuclear spins (right). 



accurate for highly polarized and moderately inhomoge- 
neous systems |25| and allows to bring results and in- 
tuition from quantum optics to bear in the spin system 
discussed here. Dark states are included in the form of 
the vacuum of the collective mode b = J^i di'^i coupled 
to the electron in Eq. ^ . For unpolarized systems (with 
on average 1/2 excitations per bosonic mode a.i), this de- 
scription provides a lower bound on the performance of 
the cooling protocol, since in the absence of an inhomoge- 
neous Knight field cooling is limited to 0(1) excitations 
per mode rather than the 0{-\fN) coolable excitations ex- 
pected at the beginning of the cooling process for spins, 
cf. Eq. (|3]). In the two limiting cases discussed so far, 
Eq. (O simplifies to 

^lij _ I -"^KSijgf^i (i) Spin Temp. 

{^ijlij - fiJ2k dkigilkj + Qolik) (m) Bosonic. 

One can take into account more aspects of the spin 
algebra by replacing some higher order expectation val- 
ues by lower orders using the properties of Pauli matrices 
[cr+,(Tj^] = o-f and afaf = icrf , obtaining 



At 



■ ^tjji] - KSij^gkigzjkj + gjjik) 

k 

k{1 - S^J)(^-J29kg^{crtcr^cr-) + 

k^i 

- 9kg J i'^t'^j'^k) + ghv 



gh^i 



(6) 



The remaining higher order expectation values (now hav- 
ing distinct indices i ^ j, j ^ k) can be approximated in 
a Hartree-like way [42] (iii), or, having the bosonic limit 
in mind, by the Wick theorem (iv), 



[jn - \)lkj [Hi), 

-\lk3 + Iktlij + Ikjlu {iv). 



The fifth and final approximation scheme we invoke 
has been introduced in the context of superradiance as 



a Wick-type factorization, that takes into account the 
partly bosonic, partly fermionic properties of spin-1/2 
operators [IHl- In contrast to the last two factorization 
schemes, it does not rely on distinction of cases. It is di- 
rectly based on the exact Eq.(l5]), and approximates the 
three-operator-expectation values in the following way 
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Ik] - Ikilij + Ikjlii (w) "Spin" 



Direct comparison of the approximation schemes (i)- 
(w) with the exact solution for both homogeneous and 
inhomogeneous couplings is shown in Fig. [21 In the ho- 
mogeneous case the spin temperature description (i) is 
clearly qualitatively wrong, because it neglects correla- 
tions in the bath. The bosonic description [ii) captures 
the feature of dark states, but it overestimates their 
influence: Instead of ~ ^fN , only one excitation can 
be removed. The two schemes based on distinction of 
cases, (iii) and (iv), give very good results initially, un- 
til roughly ^/N spins have been flipped. Then however, 
the polarization keeps increasing on a slow timescale and 
does not reach a steady state in the correct time. The (v)- 
"spin" -approximation gives very good results, and gets 
both the polarization timescale and the flnally obtained 
value of the polarization right within a few percent. 

The comparison of the different approaches to the ex- 
act solution for inhomogeneous couplings is restricted to 
small particle numbers (see Fig. [2]). In this regime all 
introduced approximation schemes reproduce the exact 
dynamics correctly. The reason for the good correspon- 
dence is the strong dephasing of dark states and gener- 
ally coherences between nuclear spins for small inhomo- 
geneous systems. 

Using these approximations we present the polariza- 
tion dynamics for N = 10"^ spins coupled through a 2D 
Gaussian wave function in Fig. [31 For the data presented 
in this and the following flgure, we considered the spins 
in a 2D square lattice geometry, with the lattice constant 
set to unity. The bosonic description displays the lowest 
final polarization and polarization rate (for the same rea- 
sons as in the homogeneous case) and is expected to give 
lower bounds on the performance on the polarization pro- 
cedure. Of particular interest are the predictions of the 
(w)- "spin" -approximation scheme, because its good per- 
formance in the completely homogeneous situation gives 
confidence that also partial homogeneities are correctly 
accounted for. Achieved polarizations of ~ 60% in this 
setting show the importance of the intrinsic dephasing 
due to the inhomogeneity (homogeneous coupling would 
allow for < 5% polarization). However, the intrinsic in- 
homogeneity alone does not allow for ultra-high polar- 
izations and we are thus lead to investigate more sophis- 
ticated cooling schemes. As shown later, in these en- 
hanced protocols all approximation schemes lead to the 
same conclusions. 

To gain a better understanding of the presented phe- 
nomena in the inhomogeneous situation, we go to an in- 
teraction picture pi = UqpUq, with Uq = exp(—iA^t/2), 
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FIG. 3: (Color online) The polarization dynamics for — 
1000 spins coupled with a 2D Gaussian wave function, which 
is shifted from the origin by 1/3 in x- and y-direction. 



which shows very clearly the oscillating coherences be- 
tween spins with gi ^ gj 



Ap 



^g^g^e~^9i9.-9:)t/2^+^j^p^ 



In the rotating wave approximation (RWA), the rotat- 
ing terms {gi ^ gj) are neglected and in the absence of 
exact symmetries the above equation reduces to the spin 
temperature description. A partial rotating wave approx- 
imation neglects only the coherences between spins with 
considerably different coupling constants, i.e. the ratio 
between dephasing and polarization rate is required to 
be large {4:\gi — gj\/ {gAtgigj) > 1). This procedure gives 
a block diagonal Liouvillian which allows for the exten- 
sion of the numerical studies to particle numbers up to 
= 10-*. 

In the RWA we evaluate the build-up time Tp for the 
polarization as the inverse of the weighted average of the 
individual spin decay times 



E.9^ J 9{9At)j:.g! \A{gAt)J 



(8) 



and find good agreement with the numerically obtained 
timescale to reach the steady state in all discussed 
schemes. For example, for the data presented in Fig. [3] we 
find times of 3.4 x 10^ (Spin Temp.), 4.6 x 10^ (Bosonic), 
and 3.3x10^ ("Spin") in units ofA''^ to reach (l-e^^) « 
0.63 of the quasi steady state Overhauser-field. This 
agrees well with the analytical estimate Tp w 2.4 x 10^ /A; 
despite the differences in the final polarizations obtained 
in the different approximation schemes. This correspon- 
dence between the RWA-based estimate and the numeri- 
cally obtained polarization times for the coherent evolu- 
tion indicates that the inhomogeneous Knight field pro- 
vides coolable excitations at a rate larger than the polar- 
ization rate, thus not slowing down the process. 



When the inhomogeneity of the coupling is large 
enough to justify the rotating wave approximation, each 
spin evolves with its own Liouvillian and the nuclei re- 
main in a product state during the whole evolution. To 
keep the errors in the derivation of the master equation 
(due to higher order terms of the expansion of the time 
evolution operators in Eq.Q) small, it is sufficient to do 
so for each spin individually in this case. This allows a 
larger time step At <C (Aamax)^^ — 0{N/A) in each 
cycle and therefore the cooling rate can be significantly 
enhanced. The cooling time effectively scales only lin- 
early in the particle number 



O 
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{jA{A/NAt)^ 



(9) 



Taking A = 100/xeV~ 40ps, a value typical for GaAs 
QDs, and 0.1 as the value for the terms gAt and A/NAt 
in the denominators of Eqs.® and ([9]) respectively, we 
find that approximately 4 x 10'^ and 3 x 10^ spins can be 
cooled to more than 90% of the steady state value {A^)ss 
within a millisecond. 

We now study enhanced cooling protocols that lift the 
dark-state limitations and which rely solely on the ability 
to shift the center of the electron wave function. These 
shifts can be effected by applying dc gate voltages to 
the QD. After such a shift only very few spins will have 
the same coupling constants for both wave functions and 
therefore singlet-like coherences are broken up. We con- 
firm this expectation numerically as shown in Fig. [4] for 
some exemplarily chosen shifts of the electron wave func- 
tion. The shifts range from a few lattice sites to roughly 
the width of the electron wave function. The timing of 
the shifts we have performed for obtaining the data pre- 
sented in Fig. m can be inferred from the plots, as it is 
accompanied by a rapid increase in the cooling rate. 

Regarding the approximation schemes, we have found 
that all schemes taking into account coherences, {ii)-{v), 
predict the same behavior, and the spin-based factoriza- 
tion (v) offers the quantitatively best description. It is 
important to note that all these descriptions coincide at 
the end of the cooling protocol [shown in Fig. 0] only for 
(ii) and (v)]. In particular the limiting bosonic model 
predicts the same high (> 95%) polarizations and cool- 
ing rates as the other schemes, which leads us to con- 
clude that 0(10) mode changes are sufficient to achieve 
near-ground state cooling for realistically large numbers 
of nuclei in QDs. 

Despite being a radical approximation at low polar- 
ization, the bosonic scheme (ii) captures the cooling dy- 
namics qualitatively and we remark that it can be gen- 
eralized to provide an accurate and conceptually simple 
description of the electron-nuclear spin dynamics at high 
polarizations [2^ . 

The cooling schemes we have presented are governed 
by the optimal timescale set by the hyperfine interaction 
constant A, but the schemes themselves leave room for 
optimization: The cooling rate can be tuned by choos- 
ing At adaptively during the cooling process. The mode 
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FIG. 4: (Color online) Polarization dynamics in the enhanced 
cooling protocol for = 196 (upper plots) and N = 1000 
(lower plot). In the upper plots approximation schemes {ii) 
(left) and («) (right) have been invoked, the lower plot is 
based on the bosonic model and the partial rotating wave 
approximation (see text). In all plots the different lines are 
representing cooling procedures with different numbers of 
modes changes. In the upper plots the randomly chosen 
Gaussian modes with width w = N/i are defined by the 
centers {(1/3, 1/3), (1.35, -0.81), (0.32, -0.04), (1.17, 0.79), 
(-0.13, -1.44), (0.96, -0.17), (0.35, 0.88), (1.27, 0.71)}. 
In the lower plot only two modes with centers 
{(1/3, 1/3), (—3.15, —1.5)} have been iterated. 



changes can be optimized by a careful choice of the size 
and the timing of the shifts, and through more sophisti- 
cated deformations of the electron wave function. These 
and further modifications are implementation-dependent 
and will be the topic of future work. 

In using the Hamiltonian Eq. ([1]) we have neglected a 
number of weak interactions that are present in actual 
systems and, while being much smaller that the domi- 
nant hyperfine term, may become important on the long 
time-scales required to reach high polarization. We argue 
in the following that these terms do not affect the quan- 
titative conclusions obtained. While nuclear Zeeman en- 
ergies are large enough to cause additional dephasing be- 
tween the nuclear spins, similar to the inhomogeneous 
Knight fields, this will only be effective between nuclei of 
different Zeeman energy, i.e., belonging to different nu- 
clear species. This leads to 2 to 3 mutually decohered 
subsystems (in a partial rotating wave approximation) 
each of which is described by our model. 

The nuclear dipole-dipole interaction [43,] can lead to 
both diffusion and dephasing processes, both of which are 
of minor importance as shown below. Dipolar processes 
that change are off-resonant and hence expected to 
be slow, as indicated by the nuclear spin diffusion rates 
measured, e.g., in [3] and should not significantly affect 



the polarizations reached. Resonant processes such as 
terms cx /f /| affect the cooling process only insofar as 
they can cause dephasing of dark states similar to the 
inhomogeneous Knight shift. The rate at which coolable 
excitations are provided is set by the energy difference for 
two nuclear spins in a dark pair. The interaction energy 
for two neighboring spins is about ~ 10~^/ieV [?!, hence 
a singlet of neighboring spins can dephase in ~ 100/zs 
(or slower if all surrounding spins are polarized). Even 
widely separated spins interacting with differently polar- 
ized environments dephase only up to a few ten times 
faster than this (depending on the geometry). Thus we 
see that the dipolar dephasing is considerably slower than 
that caused by the inhomogeneous Knight field and only 
if the latter becomes inefficient due to homogeneities 
(towards the end of cooling a given mode) the dipolar 
dephasing can contribute coolable excitations, but at a 
much slower rate than what can be achieved by chang- 
ing the electron wave function and the ensuing return to 
a situation of strong Knight inhomogeneity. Thus, one 
does not expect the cooling process to be affected ex- 
cept for a slight additional dephasing. However, on much 
longer timescales of 10s of ms the dipole-dipole interac- 
tion provides depolarizing mechanism (affecting mainly 
nuclei with a weak hyperfine interaction) that needs to 
be considered, e.g., when cooling much beyond 90% po- 
larization is studied. 

Clearly a polarization < 100% of the electron "reser- 
voir" directly translates into limitations on the final po- 
larization of the nuclei. A quantification of this neces- 
sarily needs to refer to the details a concrete physical 
realization of our model, which is not the topic of this 
article. The limitations can be minute, e.g. in the case 
of the double dot setup presented in the next section. 



IV. ADAPTING THE MODEL TO CONCRETE 
PHYSICAL SETTINGS 

The generic model of a single spin- 1/2 particle coupled 
inhomogeneously to an ensemble of N nuclear spins can 
readily be adapted to various experimental settings. 

If a source of spin polarized electrons is available, sin- 
gle electron tunneling into the QD provides the initializa- 
tion. Controlled tunneling into and out of the QD with 
rates > 10ns~^ appears feasible [H, S^, justifying the 
description of the dynamics by a suddenly switched on 
and off interaction. 

For self- assembled QDs, optical pumping with polar- 
ized light has been shown to provide a spin polarized 
bath of electrons that cools the nuclei [26l. [27l[28l[29l[47j . 
However, in this setup the average dwell time of a single 
polarized electron in the dot is large and the detuning 
due to the z-component of the Overhauser field leads to 
instabilities [H, [s^, in the nuclear polarization which 
are avoided in our scheme. 

In double QDs in the two-electron regime [1^, Hoj the 
role of the states | i), | T) is played by the two-electron 
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singlet \S) and one of the triplet states; in the following 
we consider |r+) — |t)|T)- Tunnel couphng between the 
two dots and the external magnetic field are chosen such 
that the other triplet states are off-resonant and cause 
only small corrections to the dynamics sketched here. 

As discussed in more detail in [i^ [HO, [EH the hyperfine 
interaction in this system is described by the Hamilonian 
J2i S; - A;, where I — L,R refers to the orbital state of 
the electron. Couphng between |5) and |T+) is mediated 
by the difference 6A^ = {A^ — A^)/2 of the collective 
nuclear spin operators of the two dots L, R, while the ef- 
fective Overhauser field is given by the sum {Af^ + A^j^)/2. 
Thus we have that the analysis of the previous sections 
applies to the double dot case in this regime (to zeroth 
order, cf. \E^) with the replacements 

|T)-|r+), 

A-^-V2icosO)SA-, A^^\iAl^A^^. 

The adiabatic singlet has contributions from both the de- 
localized (1, 1) and the localized (0, 2) charge states, and 
with cos9 we denote the amplitude of the (1, 1) contribu- 
tion [iOl (with (m, n) we denote a state with m electrons 
on the left and n electrons on the right dot). The effect of 
higher-order terms (e.g., of the nuclear spin components 
5A^ + A"^) merits more detailed analysis. 

This system is of particular interest since fast electri- 
cal control of gate voltages can provide a highly spin 
polarized electron system through near unit fidelity ini- 
tialization of a singlet in the right hand dot |5(0, 2)) 
[isl Issj. Starting from this singlet, rapid adiabatic pas- 
sage (Ins [ll]) by means of tuning the asymmetry pa- 
rameter e between the dots, initializes the electrons to 
the adiabatic singlet \S) and brings the system to the 
S — resonance. 

The transitions from the singlet to the other two 
triplets To._ are detuned by an external magnetic field 
(of order 100 mT in the experiments of Ref. [13|). After 
a time At the system is ramped back to the (0,2) charge 
region and the electrons relax to the singlet ground state, 
completing one cooling cycle. If relaxation to the state 
5(0, 2) is fast, the limiting timescale for this cycle is given 
by the hyperfine coupling constant A, showing that here 
the polarization rate is governed by the natural and opti- 
mal timescale (and not other, slower timescales, like e.g. 
cotunneling in Refs. [Sll. l36jl. 

In the GaAs double dot setup the sudden approxima- 
tion is justified for typical tunnel couplings ~ lO/icV, 
which have to be compared to the typical timescale for 
a hyperfine flip < O.l/icV and the fact that addition- 
ally all spin flip transitions are off-resonant during the 
adiabatic ramp. At the S — T+ resonance selecting a 
suitable combination of external magnetic field and time 
step At detunes the unwanted transitions and at the 
same time ensures resonance for the polarizing transi- 
tion. Note also that the Overhauser field increases the 
external magnetic field in materials with negative elec- 
tron (7-factor, like GaAs {g* w —0.44), thus further sup- 



pressing unwanted transitions and requiring retuning of 
the end-point of the adiabatic ramp. Given the availabil- 
ity of fast (100 ps) voltage pulses, the reinitialization of 
|S'(G, 2)) via a (0, 1) charge state is likely to be limited 
by the tunneling rate from the reservoir to the QD. For 
optimal cooling efficiency this rate should and could be 
made large > IQA/^/N |4i,|4i|. 

Since in the double dot setup the "polarized" state is 
a spin singlet, there is no inhomogeneous Knight field to 
dephase the dark states and DNP will be severely lim- 
ited. However there are many ways of providing it, for 
example by extending the cooling cycle to include a third 
step in which a single-electron state of the double dot is 
realized or by increasing the time spent at the S — 
resonance in each cooling cycle (the latter would require 
a reformulation of the master equation ([3]) not presented 
here). At the same time it would be interesting to find 
evidence for quantum coherence between nuclear spins 
in QDs by comparison of the obtained Overhauser field 
in the case of strong and weak inhomogeneous Knight 
fields [H. 



V. CONCLUSIONS AND OUTLOOK 

In summary we have presented a quantum treatment 
of a dynamical nuclear spin polarization scheme in single- 
electron quantum dots that takes into account quantum 
coherences between nuclei and allows numerical study of 
the cooling dynamics for thousands of spins. We have 
quantified limitations due to dark states and shown that 
these limits are overcome by the inhomogeneous Knight 
shift and active mode changes. From this we conclude 
that cooling to more than 90% (of the maximal Over- 
hauser field) is feasible faster than typical nuclear spin 
diffusion processes. Setups for the experimental realiza- 
tion of our scheme have been proposed. 

In order to go beyond the presented results to polar- 
izations larger than 99%, which would bring the system 
of coupled nuclei close to a pure state and significantly 
reduce electron spin decoherence, the presented scheme 
can be optimized, both in terms of timing (length of the 
individual cooling step and wave function changes) and 
in terms of the electron wave functions chosen. A further 
enhancement may be achieved by combining the polariza- 
tion scheme with -measurements [s^, IssTISGj to reduce 
the A^ variance and to tailor the interaction times and 
the external field to the measured A^ value. Dipolar in- 
teraction and other depolarizing processes will become 
more important in later stages of the cooling and need 
to be considered carefully in the development of ground- 
state cooling techniques. More detailed studies of these 
processes may, in addition, lead to schemes to monitor 
the intrinsic (dipolar) nuclear dynamics via the hyperfine 
interaction. 

The combination of high polarization and long coher- 
ence times make the nuclear spin ensemble itself a candi- 
date for an active role in quantum computation. Like the 
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actively explored single- nucleus-spin qubits [H*] , collective 
excitations of a polarized ensemble of spins could also be 
used for quantum information purposes psj . Similar to 
their atomic counterparts [stI [58| , the ensembles might 
become more suited than their isolated constituents for 
certain quantum information tasks. 
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